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g^: Abstract 



Inclusive hadro production in e annihilation processes is examined to study the fragmen- 



tation process. A broken SU(3) model is used to determine the quark and gluon fragmentation 



functions of octet vector mesons, p and K*, in a simple way with an SU(3) breaking parameter 
A. These are expressed in terms of just two light quark fragmentation functions, V^XjQ"^) and 



' 'y{x,Q'^) and the gluon fragmentation function Dg{x,Q'^). These functions are parameterized at 

the low input scale of Qq = 1.5 GeV^, evolved through LO DGLAP evolution including charm 
and bottom flavour at appropriate thresholds, and fitted by comparison with data at the Z-pole. 
The model is extended with the introduction of a few additional parameters to include a study of 
singlet-octet mixing and hence to and (j) fragmentation. The model gives good fits to the available 
data for x > 0.01, where x is the scaled energy of the hadron. The model is then applied success- 
fully to uj, (j) production in pp collisions at the Relativistic Heavy Ion Collider, RHIC; these data 
form an important base-line for the study of Quark Gluon Plasma in heavy nucleus collisions at 
RHIC, and also in future at the LHC. 

PACS numbers: 13.60.Le, 13.60.Hb, 13.66.Bc, 13.85.Ni 
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I. INTRODUCTION 



Studies of meson fragmentation are currently limited by data in comparison with the 
relatively abundant data available on deep inelastic processes (DIS). Meson fragmentation 
processes can be understood within QCD through time-hke conjugates of the space-like pro- 
cesses that contribute to DIS. Hence there is a great deal of interest in the study of meson 
(and baryon) fragmentation. Data from e+ collisions are most commonly available for 
pseudo-scalar and vector meson fragmentation. While ep data is severely limited, there has 
recently been high quality data on pseudo-scalar meson production from pp collision pro- 
cesses in RHIC[1] . Prehminary data on hght mesons is already available from the LHC[2] and 
more is expected shortly. In particular, an understanding of 77 and meson fragmentation 
m pp processes is important as a basehne for the study of the production of these mesons 
in nucleus-nucleus collisions as a signal of quark-gluon plasma (QGP)[3]. While there exist 
many phenomenological studies on tt and K meson[4-10] and 77 fragmentation[5, 11], as well 
as many comprehensive reviews[6-8, 17] of these, the issue of vector meson hadro-production 
has not been addressed so far. In this paper we focus attention for the first time on light 
{u,d, s valence quarks only) vector meson fragmentation using a model that has earlier been 
applied to a study of light pseudo-scalar meson[4, 5] and octet baryon[4] fragmentation. In 
particular, the (j) meson, which is almost a pure strange quark-anti-quark bound state, has 
special relevance as a signal for QGP [12-14]. 

A study of fragmentation functions requires an experimental input at a given (mo- 
mentum transfer) scale since QCD cannot predict the fragmentation functions themselves 
but only their Q"^ dependences. In the case of mesons formed from the light quarks {u, d, 
s), it is possible to apply symmetry arguments to reduce further the number of (unknown) 
input starting fragmentation functions. 

A simple SU(3) symmetric model is introduced[4] which has been applied to pseudo-scalar 
octet mesons and octet baryons at the leading order level[4, 5]. In this paper, the above 
model is applied to study the fragmentation functions of octet vector mesons (p and K*) 
in e"*" e~ annihilation. The model reduces the various required input quark fragmentation 
functions to a combination of just three independent fragmentation functions, a{x,Q^), 
/3(a;,Q^), j{x,Q^) and an SU(3) breaking scale-independent parameter A at a low input 
scale. In addition to these functions, the gluon and heavy quark contributions are also taken 
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into consideration during evolution. The model is then extended to predict the octet-singlet 
mixing oiuj,(l) mesons with a very simple ansatz that relates the single fragmentation function 
in the singlet sector to an octet fragmentation function along with a few additional constants. 

In contrast, other studies[6-ll] in, say, the pseudo-scalar meson sector, fit individual data 
on different mesons with no attempt made to combine the data in any way whatsoever. 
This is in fact one possible reason why the vector meson sparse data set may have not been 
studied so far. 

The paper is organized as follows: the kinematics of the relevant e"^ e~ and pp scattering 
processes is briefly explained in section 2. The model to study the fragmentation functions 
of vector mesons is introduced in section 3. The model is used to study the cross-section 
for hadro-production in section 4. In particular, a detailed study of the pure octet mesons, 
p and K*, is given in this section. It is then extended to include the singlet case, with 
singlet-octet mixing, in section 5. The detailed parameterization is explained in section 6. 
Finally, the model is compared with both e'^ e" and pp data in section 7. A summary and 
discussion is presented in section 8 of the paper. 

II. FORMALISM 

A. Hadron production in 6+ e~ process 

The e+ annihilation process is used to analyze the fragmentation functions of quarks 
(anti-quarks) through the hadrons which they produce. The reaction e'^ e~ ^ h-\- X pro- 
ceeds by flrst creating a quark and an anti-quark pair through e"^ e" — > qq via an intermediate 
vector boson V — 7/^°, followed by the fragmentation of the quark(g), anti-quark (^), or 
gluon (^f) into a hadron. This process is known as the fragmentation process[15]. 




FIG. 1: A schematic diagram of the fragmentation process in e+ e scattering[16]. 
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The term h in Fig. 1 represents the hadron (vector meson in our case), where X is the 
debris. The corresponding cross-section for this scattering process at cm. energy y/s fac- 
torises into a component describing the hard scattering and one describing the hadronisation. 

Hence, it can be expressed to LO as[17]: 

in which the fragmentation function Dg{x, Q^) is the probabihty for a quark to hadronise to a 
hadron carrying a fraction x of the energy from its parent quark, where x = Ehadron/ Equark — 
{2Eh/-\/s) < 1 (or Xp = 2ph/-\/s) and Q —y/s is the energy scale, where we are considering 
the reaction. 

The charge factors Cq are associated with the quark qi with flavor i, written[17] in terms 
of the electromagnetic charge e^, vector and axial vector electroweak couplings, Vi — T^i — 
2ei sin^ and Oj = T^i, as 

' q ' 

Cq = ^^[4 + '^^g'^eVq Pl{s) + {vl + al)v^g P2{s)] , 

s 

Cq = [v^ +aj p2{s) , 

s 

^ /X ^ 1 s{ml - s) 

' 4sin2^^cos2^w (m|-s)2 + m|r| ' ^ ' 

1 

P2[S) = 



(4 sin^ 0„ cos2 ^„)2 (m| - s)^ + m\Y\ ' 
The values of Tsj, the third component of weak isospin, and Vi are tabulated with weak 
mixing angle in [18]. Vz and mz are the decay width and mass of the Z- intermediate 
gauge boson for high energy scale. We re-express the LO cross section in terms of singlet 
and non-singlet fragmentation functions, as 

1 dcr'^ _ T,j(ijDj{x,t) 



(Ttot dx Y^qCq ' 

where j = 0, 3, 8, 15, 24 and Do, D3, Dg, D15 and D24 refer to the singlet {Dq = ^^(-C'j + D-)), 
and the non-singlet (n^ — 1 = 3,8, 15,24: {u — d), {u + d — 2s), {u + d + s — 3c) and 
{u + d + s + c — 4b); u = D„ -|- Du, etc.) combinations respectively. Here the coefficients are 
do = (cu+Cd+Cs+Cc+Cb) /5, 03 = {cu-Cd)/2, as = {cu+Cd-2cs)/6, 015 = {cu+Cd+Cs-3cc) / 12 
and = {cu + Cd + Cg + Cc — Acb) /20. Note that the data from LEP or SLD on the Z- 
pole is dominated by the Dq fragmentation function combination while at lower energies 



where photon exchange dominates, the u-type quark fragmentation functions dominate. 
Furthermore, data is available[26-29] only for the sum of the conjugate states, — p'^ + p~, 
and similarly for the charged and neutral K* mesons [3 0-3 3]. We therefore compare with 
the sums of these quantities. The SLD data[40, 41] is available separately for fragmentation 
from light quarks (m, d, s) alone; hence for these data, the summation in Eq. (1) is over the 
light flavour contributions only. Data on the singlet-octet mixture of states, viz., ou and (f), 
are available from both LEP and SLD [34-41]. 



B. Hadron production in pp process 

In addition to information on the unknown fragmentation functions, hadro-production 
in pp processes requires a knowledge of the parton (both quark and gluon) distribution 
functions within the primary proton. The scattering is expressed in terms of the underlying 
parton interactions, with one of the final state partons fragmenting into the hadron of 
interest. All possible initial state parton interactions are possible: qq, qg and g g, as well 
as processes with antiquarks. 

The inclusive cross section for hadro-production, p + p ^ h + X , producing a hadron h 
at large pt is given by[19-21], 

where the sum over (a, b, c, d) runs over both quarks and gluons. Here Xa and Xb are the 
usual Bjorken-x variables corresponding to the parent proton momenta p^ and ps'- Xa = 
Pa/PA, = Pb/PB (neglecting intrinsic transverse momentum) and Pa/A{xa,Q'^) are the 
usual parton density distributions; for example, Pu/p{xa, Q^) = u{xa, Q^), etc. 

The fragmentation functions depend on the variables, z — Zh — Ph/Pc, the fraction of 
momentum of the quark carried by the fragmenting hadron, and the scale ^ Pt- The 
limits of integration are[19-21] 

" 1-X2 ' Xa-Xi ' 

with xi — —u/s,X2 — —t/s. 

For numerical comparison with the data, we reexpress the cross-section in terms of the 
physical observables which are the transverse momentum pt — Ph sin 9 and the rapidity 



y^{l/2)\n[{Eh+PhC0se)/{Eh-phC0se)], as 

^ dV _ 1 dV 
^dpl ~ prdprdydcl) ' 

where 6 is the scattering angle of the hadron h in the pp center of mass frame and Eh and 
Ph are its energy and 3-momentum. Note that the azimuthal angle dependence is trivial in 
this process; furthermore the data from the PHENIX experiment at RHIC is taken over a 
range tt in and —0.35 < y < 0.35 in the rapidity [12-1 4] . The sub-process cross-sections 
are well-known[19-22]; the qq, qg and gg processes all contribute at the same order in ccg. 
Hence the quark and gluon fragmentation functions contribute at the same order, unlike 
in the e+ e~ case. This data is therefore an important test of the correctness of the gluon 
fragmentation functions. 

We now present details of our model for quark fragmentation functions. 



III. SU(3) MODEL 

Fragmentation functions parameterize the hadronisation process, in which the observed 
hadrons are formed from the final state partons of the scattering process; these cannot 
be calculated in QCD. However, given a definite energy scale and starting distribution, 
QCD can evolve these perturbatively and explain their scale (Q^) dependence. A set of 
common fragmentation functions are used to describe the members of the octet of vector 
mesons K* {K*+ , K*' , K*'^ , K ) and u. An SU(3) symmetric model (with 

parameterized SU(3) breaking) has been chosen to achieve this because, in principle, such 
a symmetry gives good description about the octet of vector mesons. Such a model with 
broken SU(3) was developed for pseudo-scalar octet tt and K meson fragmentation in e"*" e~ 
collisions[5] and for the octet baryons[4] p, n, A, E. We use an analogous model in the vector 
meson sector in our analysis. 

We start with light quarks {u, d and s) at the starting scale of Qq — I.^ GeV^. However, 
as the fragmentation function evolves across various thresholds (typically up to = (91.2 
GeV)^), consistent contribution of charm and bottom quarks are included appropriately in 
the evolution. Let us consider the process at the input scale as, 

qi h] + Xj . 
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The underlying SU(3) process can be thought of as 3 — > 8 + X. That is, a quark goes 
to an octet hadron (/i*) with the remainder Xj being a triplet (3), antisixplet (6) or fif- 
teenplet (15) with i,j running over 1 to 8. Let Q;(a;,Q^), /3{x,Q'^) and 7(x,Q^) be the 
corresponding unknown SU(3) symmetric independent fragmentation functions for each of 
these possibilities[4], that is, for X to be 3 (6, 15) the probability of the quark to fragment 
into an octet meson is a {f3,j). 

In a similar way, an anti-quark also produces an octet hadron with Xj being an anti- 
triplet (3), sixplet (6) or anti-fifteenplet (15), for which a{x,Q'^), /3(x, Q^) and 7(^,(5^) 
have to be determined. Thus a single meson has seven unknown fragmentation functions 
D^{x,Q'^), D^-^{x,Q'^) and Dg{x,Q'^) associated with its production. Here D'^, and 
refer to the light quark, anti-quark and gluon fragmentation functions while the heavier 
quark contributions are zero at the starting scale (below the charm threshold). 

So, we have to fit a total of 56 (8 x 7) unknown fragmentation functions to the data for 
octet mesons, which is rather daunting. The problem is made simpler when we apply SU(3) 
symmetry since the seven independent symmetric fragmentation functions a{x, Q^), f3{x, Q^) 
and ■j{x,Q^) including their conjugates a(x, Q^), f3{x,Q'^) and 7(0;, Q^) and Dg[x,Q'^), the 
gluon contribution, should determine the fragmentation of the entire group of octet mesons. 
Since SU(3) symmetry is only approximate, we use a single x- independent parameter to sig- 
nify SU(3) breaking, while SU(2) remains unbroken in our model. Hence the fragmentation 
functions of and p° (and similarly for isospin conjugates of K*) are related by isospin 
symmetry. 

A. Valence and Sea functions 

The vector meson octet is a self conjugate octet. So, = . Therefore, we have three 
independent fragmentation functions as mentioned in the above section and we express 
the quark fragmentation in terms of these three fragmentation functions in Table I. We 
reduce the number of unknown functions further through various symmetry considerations 
like isospin invariance and charge conjugation. We assume the sea is flavour symmetric, so 
that DP~ = and so on. Using this assumption and the expressions for corresponding 
fragmentation functions given in Table I, we have 

^{x,Q') = ^{x,Q')/2; (7) 
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and all the sea fragmentation functions are equal to 

5(x,Q2) = 27(x,Q2) . (8) 

Thus, all valence fragmentation functions can be expressed in terms of the function V{x, Q'^), 
where V is given, say for p+, by the difference D^"*" — D^. Therefore, we have 

Vix, Q^)=(a + I3 + ^7) Q') - 27(2^, Q') , (9) 

Substituting the value of (3 from Eq. (7) in Eq. (9) we get the valence combination to be 

V{x,Q') = a{x,Q')-^^{x,Q') . (10) 

All sea fragmentation functions can be expressed in terms of 7. Thus, Eqs. (8) and (10) 
represent just two unknown fragmentation functions for the sea S (or equivalently 7) and 

valence V, in terms of which all quark fragmentation functions which describe all the octet 
vector meson production can be represented. 

B. Breaking of SU(3) symmetry 

SU(3) symmetry is broken in the model due to relatively more massive strange quarks. So, 
in addition to valence and sea functions, we introduce an x-independent symmetry breaking 
parameter A for a non-strange quark to fragment into a strange octet meson. 

For example, to produce K*~^ (us) meson, a u quark in the valence has to pick up a(n 
anti) s quark; being more massive, the corresponding fragmentation function is suppressed 
by the parameter A. Moreover, if the valence quark is s then it may easily pick up a w quark 
without this suppression factor. In the same way, fragmentation of other mesons like K*~ , 
K*'^ and K*^ can be explained. 

For all these mesons with strange quark in their valence, apart from different valence 
parts, their sea function remains the same — 2A7 — because the sea is flavour symmetric and 
is uniformly suppressed by the suppression factor A. Thus we use broken SU(3) symmetry 
with an SU(3) symmetric sea as our model to describe the fragmentation functions of vector 
mesons with the introduction of the parameter A. 
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IV. FRAGMENTATION OF p AND K* MESONS 



We begin by neglecting the problematical uj meson that is not a pure SU(3) octet meson; 
it will be considered along with the meson later. 



A. p meson 



The scattering cross section is expressed in terms of singlet and non singlet combinations 
in Eq. (3). Let us begin with p, the lightest meson, and express the fragmentation function 
for u quark from Table I. On substituting the values of sea and valence parts from Eqs. (8) 
and (10) we get, 

L'f =y + 27, 

as u quark is present both in the valence as well as in the sea part and is equal to the d 
contribution. The charge conjugation invariance clearly shows that the other quarks are 
not in the valence of the p meson; hence their fragmentation for sea part has to be 27. 
The singlet contribution at the input scale for p meson (where only the three light flavours 
contribute) is therefore the total quark contribution: 

Dg+ = D''*_ - = 2V + 127 . (11) 

U u+u+d+d+s+s ' ^ ' 

The non singlet contributions D3 and can be obtained for p meson with the same analogy. 
Here, the non singlet term D3 (= — turns out to be zero due to charge conjugation 
invariance and — 2V. 



B. K* meson 



The procedure for the K* meson is exactly the same as above. The only difference is that 
it has a strangeness quantum number. Thus for the valence component of K*" meson (su), 
the non-strange u quark contribution is suppressed by A as it has to pick up a strange quark 
to form the K*~ meson. However, the s contribution is itself not suppressed, since only an 
(anti-)M quark is required here. Meanwhile, all the sea fragmentation functions come with 
a uniform factor of A since for all quarks (flavour symmetric sea) , a strange quark has to be 
produced. Therefore, 

Dr - D-:^i^d.-d^s^-s = (1 + >^)y + 12A7 . (12) 
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Notice that while A is x- independent, there is an inherent x-dependence of the strangeness 
suppression, being maximal at small-x and least at large-x due to the different suppression 
of the valence and the sea quarks. 

Thus p and K* mesons have different valence and sea quark fragmentation functions. In 
addition, since the gluon fragmentation function mixes with the singlet Dq fragmentation 
function on evolution, we parameterize a possible gluon suppression through D^* — f^* D^. 

V. EXTENSION TO u AND </> MESONS 

We now extend this broken SU(3) model which explains in a most simple way the pure 
octet (p and K*) mesons, to the ui and (f) mesons, which are orthogonal combinations of the 
SU(3) octet (cjg) and singlet states (cui): 

ui — sin ^ a;8 + cos 6 uji , 

(f) — cos9 (jj^ — sn\6 uji . (13) 

where = {uu + dd - 2ss) /Te, ooi = {uu + dd + ss) I \/3 are the corresponding orthogonal 
states and d is the vector mixing angle, whose value is approximately[23] 35°. Note that 
a value of Q close to this value saturates the physical ^ state[24] as a pure ss state. The 
fragmentation functions for uj^ can be described with the help of fragmentation functions 
given in Table I, since uj^ is one of the members of the octet, whereas details regarding o^i 
will be discussed in the next section. 

A. Singlet hadron (cji) fragmentation 

Let us consider the same process, 

in which a quark hadronises into a singlet meson so that X can only be a triplet (3 — )■ l+X). 
Therefore, we need to determine only one unknown fragmentation function b[x^ Q"^) in the 
singlet case. In section 3 we saw that the probability for a parton to fragment into an octet 
hadron with X being triplet is a{x, Q"^). For fragmenting to a meson so that the remainder 
X is a triplet state, therefore, there are only two possibilities: either the hadron is an octet 
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(the process is proportional to a) or the hadron is a singlet (the process is proportional to 
5). Hence we use the simple ansatz that the function 5 is simply related to a{x,Q'^), the 
fragmentation function for members of octet meson. That is, 



where the factor of 1/3 is due to the normahzation of the state and /i is the (presumed) 
a;-indcpcndcnt proportionality constant we have to determine in the analysis. Since this 
uses an approximate SU(3) symmetry, this constraint is applied only at the input scale Ql 
where there are only three active flavours. The evolution is correctly applied to all the active 
flavours, depending on the scale. 

We will now express the fragmentation of uj and (f) mesons in terms of the SU(3) octet 
and singlet contributions. 

B. Singlet and octet contribution to oj and (p mesons 

The strangeness suppression factor A remains the same here. However, there is a difference 
in the sea suppression factor. While it was also equal to A in the case of the K* meson, 
here the suppression is rather different. It arises due to the preference of a given qq pair 
to fragment to the lighter p rather than to a; or 0. In particular, in the case of ou, which 
is mostly saturated by non-strange quarks in the valence, the sea suppression factor arises 
due to the preference of a given light qq pair to fragment to the lighter p rather than the 
UJ. Hence the sea suppression factor is expected to be f^^^ ~ m'^/m'^ just as in the case 
of pseudoscalar mesons[5], and we do not expect a large suppression. On the other hand, 
the physical state is saturated by the strange contribution; hence the suppression factor 
/s^a, 4> is expected to be close to A^ as a quark has to pick up both s and s from flavour 
symmetric sea. We will see in the next section that numerical fits to the data do indeed 
agree with these expectations. 

We first explicitly write down the fragmentation functions given in Table I for corre- 
sponding octet hadrons in terms of valence and sea sectors, including the various suppression 
factors as discussed above: 




(14) 




(15) 




-XV + 2/sea7 ; 
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where the sea 7 and valence V have their usual definitions, described in Eqs. (8) and (10), 
A is the strangeness suppression factor while /sea is the unknown suppression factor for the 
SU(3)-symmetric sea fragmentation functions. 
Using our ansatz for the singlet hadron, we have 

Dl-Dl^ll-(v + ^/,ea7) , (16) 

Dl = f + ^/sea7) . 

Here we have introduced separate suppression factors for the u, d- and s-type singlet frag- 
mentation functions. With these four equations in hand, we express the fragmentation 
functions for uj and mesons, given the definitions of the states in Eq. (13), with vector 
mixing angle 9, at the input scale as, 

D* = (c*)^ (cos^ + W 0^) ; (17) 

Here, i refers to the three light quarks {u, d, s); the co-efficients are cf^ = = {cos 9 — 
\/2sm.9), cf = (-2 cos 6* - a/2 sin 6*) and 4 = ^, 4 = 2, 4 = = a/2. Obviously, we can 
find the coefficients for ou meson in the same way [5]. These fragmentation functions can be 
re-expressed in terms of Dq{x,Q^), etc., in the usual way. Finally, we again parameterize 
the gluon fragmentation functions as Dg''^ — fg''^D^. 

VI. PARAMETERIZATION OF THE INPUT FRAGMENTATION FUNCTIONS 

The unknown functions for the valence V{x, Q'^), sea 7(2;, Q'^), and gluon Dg[x, Q^) frag- 
mentation are parameterized at low input scale of Ql — 1.5 GeV^ for three light quarks {u, 
d, s) where the charm and bottom contributions {Dc, D^) are zero. The parameters are then 
determined through fits with data. We use a standard functional form to describe these 
quantities: 

Fi{x) = aiX^'{l - xY'il + diX + Cix'^) , (18) 

where a^, bi, Ci, di and Cj are the values to be determined. The fragmentation functions 
are evolved to leading order (LO) including charm and bottom contributions in appropriate 
places along with the gluon fragmentation, to the Q'^ values of interest. Since the gluon 
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fragmentation does not occur in the expression for the cross-section, it is least constrained 
by these fits and relatively unknown. Hence we set eg — for the gluon. 

VII. DATA ANALYSIS AND RESULTS 
A. Fragmentation in e"*" e~ process 

We therefore have a set of paramctcrizations as well as a bunch of scale-independent 
constants that need to be determined. We focus here entirely on the data as the 

cleanest sample. Hence necessarily the fits to the gluons are ill-determined, since, to LO, 
the gluon contributes only through evolution. Furthermore, we concentrate on the Z-po\e 
data[25] from both LEP as well as SLAC-SLD. At this scale, five quark (qq) flavour pairs 
from utob are produced in the final state and hence contributions to the final state meson are 
of two types: first, is the direct fragmentation of the quark (anti-quark) into a vector meson; 
another is the fragmentation of a heavy (for example D or B) meson which subsequently 
decays into one of the mesons of interest. Since b quarks decay dominantly (greater than 
95% of the time) into c-quarks, the decay chain of both heavy b and c end up in strange 
mesons, driven by the large Vcs CKM matrix element. We therefore expect that data on the 
strangeness containing K* and mesons will have large contamination from these heavy 
flavour decays. 

On the other hand, the K mesons decay predominantly into tt and 77; hence there is very 
little contamination of the p and cu data. In this analysis, therefore, we consider the inclusive 
hadro-production data from LEP on the predominantly non-strange vector mesons, p and 
a; [24] to arise from hght quark fragmentation. For the strange mesons, K* and 0, we use 
instead the so-called uds data from SLD where the fragmentation of the light quarks alone 
have been segregated and studied. 

1. For pure octet mesons, p and K* 

For extraction of fragmentation functions from the data we used inclusive particle pro- 
duction in e" as the process of choice. The unknown fragmentation functions V{x,Q'^), 
7(2;, (5^), Dg{x,Q^) and the suppression factor A are determined by comparison with the 
data[25]. 

13 



The comparison with data is restricted to the range x > 0.01 since it is well-known that 
the DGLAP evolution does not explain the behaviour of the cross-section at small-x and 
modifications such as MLLA (modified leading log approximation) are needed to explain the 

small-a; data. 

Note that individual data sets are some-what small, with typically 6-8 x-bins. This is 
insufiicient to fit all the unknown fragmentation functions (quark and gluon). However, the 
model can be applied to the entire octet meson data as a whole; this allows for a reasonable fit 
to be obtained with available data. This, in fact, is the primary motivation for constructing 
such a model. 

Using the available data we fit the functional form of Eq. (18) for pure octet mesons 
(p and K*) by fixing the values of parameters a, 6, c, d and e for valence, sea and gluon 
fragmentation functions. 

The smaller x behaviour is dominated by the sea (7) contribution and large x behaviour 
by the valence V contribution; this helps determine the values of a, 6, c. This process is not 
very sensitive to the gluon contribution since this contributes only in the evolution and not 
directly in the definition of the cross-section at LO, unlike in, say, the pp process. 

a. The fragmentation functions : The best-fit values of the parameters in the input 
fragmentation functions and the l-cr errors on them are tabulated in Table II. These cor- 
respond to the input fragmentation functions for p mesons having three light quarks (the 
heavier c, b quarks do not contribute) at the starting scale of Qq = 1.5 GeV^, as shown 
in Fig. 2. This figure clearly shows that the valence contribution at large x is dominated 
by u and d fiavours, while strange fiavour and gluon contribute only in the small x sea 
part. As the evolution crosses various thresholds the charm and bottom contributions are 
included appropriately as depicted in Fig. 2 which clearly reflects the (small) charm and 
bottom contribution after evolution at the Z-pole. The heavy fiavour contribution purely 
arises from gluon-initiated processes and cannot account for the contribution in inclusive 
hadro-production data arising from production and decay of heavy fiavour mesons. 

h. The p, K* cross-sections at the Z-pole : Given these fragmentation functions, finally 
we plot the various cross-section combinations as given in Eq. (1) along with data [26-28] 
for the p meson in Fig. 3. The values of the fits (obtained by averaging the cross-section 
over each bin and comparing with the data) are tabulated in Table IV. (The difference 
between the average cross-section and its value at the average x value of the bin indicated 
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the uncertainty due to bin width and has been included as an error in the computation. 
This affects the significantly only in the last bin). 

We also use the fitted values to predict and compare with data on p meson fragmentation 
at the photon-exchange-dominated regime of ^/s = 29 GeV; see Fig. 3. Note that the 7 
exchange process at lower energies is sensitive to a different combination of the fragmentation 
functions than at the Z-pole. 

The value of the suppression factor can be determined by the fraction of p and K* meson 
at small x (the data at large x have relatively larger error bars). A fit to the data as seen in 
Fig. 4 gives A = 0.063 (see Table III). Notice that the 1-a range of A (0.05 < A < 0.07) for 
the vector meson octet is close to that obtained (A ~ 0.08) for strangeness suppression of 
the pseudo scalar mesons K^, and K'^. This may indicate that the origin of strangeness 
suppression may be independent of the spin structure of the mesons, since the two octets 
are otherwise unrelated. 

The data[40, 41] for K* meson is also plotted along with our predictions, using the same 
input parameters as used for the p meson, in Fig. 4 with the best fit value of the suppression 
factor for gluon, f^* = 1.0 (error bar is in Table III). The data are reasonably well fitted 
over a large x region. Thus all the parameters in the input fragmentation functions are 
completely determined from fits to the p and K* meson data. The p data gives a set of V , 
7 (and Dg), which consistently fit the K* data with the inclusion of the two parameters A 




FIG. 2: (L) Three input flavours u = d and s with gluon g at the starting scale of ^/s = \/L5 GeV, 
as a function of x, with zero contribution of c and b flavours for p meson. (R) Contribution of all 
the five flavours u = d, s, c and b with g after leading order evolution at a scale of y/s = 91.2 GeV 
for p meson. 



15 



and fg . Note that A contributes differently at low and high x: the suppression factors are 
(1 + A)/2 and A for the valence and sea quark fragmentation functions compared to those 
for p fragmentation, as seen by a comparison of Eqs. (11) and (12). Hence the excellent fits 
to the p and K* meson data validate our simple model in a non-trivial way. 

Thus using the available data for p and K* mesons we have fitted the input fragmentation 
functions as listed in Table II using a simple broken SU(3) model. We now go on to apply 
the model in the singlet-octet mixed u-(f) sector, with the fragmentation functions V, 7, and 
A fixed at the values obtained from fitting the p, K* data. 



2. For mesons of octet and singlet mixture 

As our model explains pure octet mesons in a nice fashion, we extended it to mesons 
having mixture of octet and singlet components with same fragmentation functions V, 7, 
Dg and breaking parameter A. An additional fragmentation function expresses the singlet 
contribution. However, this is simply related to a known octet fragmentation function, as 
discussed earlier, and hence only a few additional constant parameters are introduced in this 
extended sector, namely /f , fsea and fg for u and 0. 

These describe the suppression in the sea as well as relate the singlet fragmentation 
function to the octet one. With the help of Eqs. (16), (17) and (18) the fragmentation 




p "p 

FIG. 3: Fit for rho meson in terms of fragmentation functions with (L) LEP data on the Z-pole 
and (R) HRS data at ^ = 29 GeV. The data[26-28] at the Z-pole and from HRS[29] are shown 
with statistical and systematic errors added in quadrature while the dots in the left-side plot show 
the fit when averaged over the same x bins as the data. 
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functions for u and mesons are expressed in terms of those of the and UJi mesons, 
including these unknown parameters. The first is that of the mixing angle, 9. It is known 
that 9 is large and positive[24], close to 35°. Exactly at 35°, the u meson is almost purely 
non-strange while the (f) is almost purely a strange ss hadron. A simultaneous best fit to uj 
and (f) data (keeping the V, 7 and A fixed to the best fit values from the p and K* analysis) 
gives 9 = 42.6° (error bar is in Table III), not far from maximal. Hence the strange (non- 
strange) quark contribution to uj{(f)) is highly suppressed (the coefficients (c^ )^ and (c^)^ are 
just a few percent of (c^)^ and (cf)^ respectively.) Note the data are inconsistent with no 
mixing, 9 = 0°. 

Hence the u meson is totally dominated by u and d light quarks. So it only has contri- 
bution from these two fiavours whereas strangeness has least contribution. Therefore, we 
expect this to be similar to p, as is borne out by the similarity in cross-section behaviour. 
Since the strange component of u is highly suppressed, we fix the strange singlet suppression 
factor to be /f = for u. We find the data fit best to a sea suppression factor f^^^ = 0.94 and 
fg=1.0 with /" consistent with zero (error bars are in Table HI). This clearly shows that 
u behaves like p, with unsuppressed sea and gluon fragmentation functions and very little 
contamination to the non-strange fragmentation functions from the singlet mixing, that is, 
the singlet contribution merely serves to make u a practically non-strange meson orthogonal 
to p. 

Similarly, we set = for (f). Since (f) is mostly saturated by the strange contribution, 

102 
101 

Q. 
X 

I 10° 
10-1 

10-2 

0.01 0.1 1 

FIG. 4: Fit for K* meson with the best fit value of suppression factor A = 0.063. Data[40, 41] at 
the Z-pole, from light quarks only, are also shown. See caption of Fig. 3 for more details. 
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the sea fragmentation function involves picking up both s and s quarks and we therefore 
expect the sea suppression factor to be of the order of A^. Since the best-fit resuhs were 
close to this value, with large error bars, we simply set /^^^ = A^. Meanwhile, the best fit 
value of the singlet constant is /{* = 4.0. The gluon suppression factor for (p is = 0.32 
(error bars are in Table III). Hence, in contrast to the other member of the nonet, is 
heavily suppressed both with respect to sea quark and gluon fragmentation. These best-fit 
values result in a cross-section behaviour as shown in Fig. 5. 

A detailed analysis of the gluon contribution and its relative suppression in mesons is 
best done in a next-to-leading order (NLO) analysis where the gluon fragmentation func- 
tion directly appears in the expressions for the cross-section. An alternative is to study 
hadro-production in pp scattering, where gluon fragmentation is equally dominant as quark 
fragmentation due to the type of parton-level processes involved, as stated earlier, although 
there are extra uncertainties due to convolutions with parton density distributions in the 
cross-section formulae. We therefore apply the fits obtained from hadro-production in e"*" e~ 
processes, with a clean (non-hadronic) initial state, to hadro-production in pp collisions. 



B. Pragmentation in pp process 

The PHENIX experiment at RHIC has measured u and vector meson production[12- 
14] in pp collisions at ^/s = 200 GeV as a function of the transverse momentum, px- We 
compute the relevant hadro-production cross-section as expressed in Eqs. (4) and (6). We 




p 

FIG. 5: Fits to omega (L) and (R) phi meson. The data correspond to LEP data[34, 35] for co and 
the SLD data[40, 41] from light quarks afone for 0. See the caption of Fig. 3 for more details. 
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integrate over a range tt in the azimuthal angle and over a rapidity range —0.35 < y < 0.35 
as stated earlier, and compare the cross-sections, computed at a scale = Vr^ with the pr- 
dependent data. We use (GRV98-LO)[42], a standard set of parton distributions as available 
in the CERN-libraries; a different choice of parton distributions will not affect the results 
since the (x, Q"^) range of the data are in the well-studied range. The data are binned in 
Pt and the cross-sections are quoted at the central value of the bin. As with e"*" e~ data, 
the difference between the average cross-section and the cross-section value at the average 
Pt of the bin is included in the error while computing the of the fits. Note that the 
fragmentation functions are taken from the fits to the e~ data and there are no more free 
parameters. 

The results of the computation are shown for both u and (j) mesons in Fig. 6 in comparison 
with data for which the scale is larger than the starting scale of evolution, Pr > Q})- The 
figure also shows the band due to the scale uncertainty over the range p%/2 < Q"^ < 2p\,. 
It is seen that the model provides a good fit to the data. The values of are listed in 
Table IV corresponding to the central value, = p^. 



VIII. SUMMARY AND DISCUSSION 



Fragmentation functions of quarks and gluon for octet vector mesons are studied for 
the first time in inclusive hadro-production in electron-positron annihilation process and 




6 7 8 9 10 11 12 2 2.5 3 3.5 4 4.5 5 

(GeV/c) (GeV/c) 



FIG. 6: Cross section as a function oipT for omega (L) and phi (R) meson hadro-production in 
collisions. The data correspond to PHENIX/RHIC data[12-14] for uo and (p. Bands show the scale 
uncertainty on changing = over a range p\/2 (upper curve) <Q^< (lower curve). 
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proton-proton collisions to the leading order in QCD. 

The unknown fragmentation functions were fitted to data on the Z pole in e~ col- 
lisions. The model fairly depicts the octet mesons p, K* and ui and greatly reduces the 
number of unknown fragmentation functions since the symmetries reduce the 48 quark frag- 
mentation functions of the meson octet to combinations of just two: the valence V{x.Q'^) 
and sea 7(2;, Q^) fragmentation functions, apart from the gluon Dg{x,Q^) fragmentation 
function. This gives the model great predictability, since the various (sparse) data sets can 
be effectively combined to improve the quality of the fits. 

An SU(3) breaking parameter A was introduced for the strange K* mesons phenomeno- 
logically. to account for strangeness suppression, at the starting scale of evolution. The 
parameters at the input scale were determined by comparison with data on the Z pole. 
Heavy quark contributions that are radiatively generated, are small and seen to contribute 
at smaller x values. The input distributions were tuned to give good fits with inclusive hadro- 
production data at LEP for p mesons and the "pure uds" (tagged jet) ir*-production data at 
SLD on the Z-pole. The best-fit values, along with the 1-a errors are listed in Table II. The 
best fit value of the strangeness suppression factor is 0.05 < A < 0.07, close to the similar 
suppression factor obtained in fits to the pscudoscalar octet mesons, A ~ 0.08, indicating 
that strangeness suppression in quark fragmentation functions may be a spin-independent 
phenomenon. 

The model was extended with some further assumptions to include a study of the singlet 
vector meson with singlet-octet mixing (that is, to the ou and (f) mesons). While the mixing 
angle is known to be close to ^ ~ 35°, its value was determined through a simultaneous best 
fit to the u and data to be ^ ~ 43° (see Table III). This is ultimately the focus of this work 
since rj and hadro-production in pp collisions are used as a baseline when studying the 
production of these mesons as a signal of quark gluon plasma in nucleus-nucleus collisions. 

In order to describe the singlet sector, an ansatz was made, relating the singlet frag- 
mentation functions to one of the octet fragmentation functions {a{x,Q^)). The constants 
of proportionality, ff, and the gluon suppression factors were determined through the fits. 
While extending this model to include the octet and singlet nonet, no new fragmentation 
functions were introduced: just a few more parameters were included as described above, 
as well as the parameter fsea to describe suppression of sea quark fragmentation in both uj 
and (especially) (f). 
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The proportionality constants /" and /f and the sea suppression factor /sea were deter- 
mined by fitting the functions with the data. Values of the parameters obtained appear 
to be phenomenologically reasonable and are listed in Table III, along with the l-cr errors, 
while the of the fits are given in Table IV. 

Note that the input fragmentation functions at a low energy scale, {Ql = 1.5 GcV^) were 
included for three light fiavours (m, d, s) only. The heavier fiavours charm and bottom are 
consistently included at the appropriate thresholds during (leading order) evolution to the 
Z pole. As is well known, next-to-leading order corrections will worsen the fits at small-x 
unless mass corrections are added. This is beyond the scope of the present work. 

Note also that the contribution at the Z-polc is dominated by the flavour singlet fragmen- 
tation function Dq^ which is the sum of the individual quark flavour contributions with equal 
weight. Hence the fits are most sensitive to this sum rather than to the individual fragmen- 
tation functions. However, within the model, there are only two independent fragmentation 
functions, a valence and a sea combination, in terms of which all quark fragmentation func- 
tions are expressed. These appear with different weights in the cross-sections for p and K* 
production and hence can be individually determined. 

At LO, the fits arc rather insensitive to gluon fragmentation, although it appears that the 
gluon fragmentation of is significantly suppressed related to the others. Also, we have not 
included any isospin breaking effects. In particular, charge asymmetries in fragmentation 
functions are best studied through fragmentation in ep scattering, while pp processes are 
sensitive to the gluon fragmentation. However, data, as well as analysis, is not as clean 
in this sector due to uncertainties in scale, for instance, in the pp case. Hence fits to the 
fragmentation functions from data in sector can then be used as constraints while 

analyzing data from these other processes. Such an analysis was done for oj and hadro- 
production in pp collisions and was found to be in good agreement with RHIC/PHENIX 
data. Note that there are no free parameters in this fit since all the fragmentation functions 
are determined from earlier fits to the LEP data. Reasonable values of were obtained, 
as fisted in Table IV, although the scale dependences are quite severe. This is expected 
to improve at next-to-leading order, although again this is beyond the scope of the present 
work. 

In summary, fragmentation of the entire nonet of vector meson is explained via a simple 
model with broken SU(3). The model includes drastically few fragmentation functions and 
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some constant parameters at a low input scale which were then evolved to the scale of the 
data (mainly at the Z-pole). While the fits to the pure octet vector mesons were very 
good, reasonable fits were obtained when the model was extended to study the mixed ou 
and mesons. The model continued to give good fits to the hadro-production data in pp 
collisions as well, where the gluon fragmentation becomes important. This reflects the great 
predictability and efficiency of this model, especially in view of the paucity of data in this 
sector. 
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TABLE I: Quark fragmentation functions into members of meson octet in terms of tlie SU(3) 
functions, a, (3 and 7. 
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TABLE II: Best fit values of the parameters defining the input fragmentation functions at the 
starting scale of = 1.5 GeV^ , with their 1-a error bars. 
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TABLE III: Best fit values of the parameters defining the input fragmentation functions at the 
starting scale of = 1.5 GeV^, with their 1-a error bars. 



Central Value Error Bars 
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TABLE IV: for fits to inclusive vector meson production data from e experiments on the 
2^-pole from LEP and SLD experiments and from pp experiments from PHENIX at RHIC. 
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